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ABSTRACT 

Recent XMM-Newton studies of X-ray variability in the hard states of black hole X-ray 
binaries (BHXRBs) imply that the variability is generated in the 'standard' optically- 
thick accretion disc that is responsible for the multicolour black-body emission. The 
variability originates in the disc as mass-accretion fluctuations and propagates through 
the disc to 'light up' inner disc regions, eventually modulating the power-law emission 
that is produced relatively centrally. Both the covariance spectra and time lags that 
cover the soft band strongly support this scenario. 

Here we present a comparative spectral-timing study of XMM-Newton data from 
the BHXRB SWIFT J1753.5-0127 in a bright 2009 hard state with that from the 
significantly fainter 2006 hard state, to show for the first time the change in disc 
spectral-timing properties associated with a global increase in both the accretion rate 
and the relative contribution of the disc emission to the bolomctric luminosity. 

We show that, although there is strong evidence for intrinsic disc variability in the 
more luminous hard state, the disc variability amplitude is suppressed relative to that 
of the power-law emission, which contrasts with the behaviour at lower luminosities 
where the disc variability is slightly enhanced when compared with the power-law vari- 
ations. Furthermore, in the higher-luminosity data the disc variability below 0.6 keV 
becomes incoherent with the power-law and higher-energy disc emission at frequen- 
cies below 0.5 Hz, in contrast with the coherent variations seen in the 2006 data. We 
explain these differences and the associated complex lags in the 2009 data in terms of 
the fluctuating disc model, where the increase in accretion rate seen in 2009 leads to 
more pronounced and extended disc emission. If the variable signals are generated at 
small radii in the disc, the variability of disc emission can be naturally suppressed by 
the fraction of unmodulated disc emission arising from larger radii. Furthermore the 
drop in coherence can be produced by disc accretion fluctuations arising at larger radii 
which are viscously damped and hence unable to propagate to the inner, power-law 
emitting region. 

Key words: accretion, accretion discs — black hole physics — stars: individual 
(SWIFT J1753.5-0127) — X-rays: binaries 



1 INTRODUCTION 

Variability on a broad range of time-scales, ranging from 
milliseconds to hours, is a well-known characteristic of black 
hole X-ray binaries (BHXRBs). The short-term variability 
(time-scales up to a few hundred seconds) is strongly depen- 
dent on the spectral state of the object, being particularly 
enhanced in the so-called hard state, with rms amplitudes 
up to ~ 40 per cent of the average flux (Belloni ct al. 2005 
| Remfllard fc McClintock|2006nMunoz-Darias, Motta, fc Bel-| 
|loni|2011| |. 
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In the hard state, the X-ray emission is dominated by 
a hard power-law component (typically, T ~ 1.4 - 2.1) 
accompanied by weaker multi-colour black-body emission 
associated with the accretion disc, with fcTjn < 0.5 keV 
(Miller, Homan, & Miniutti 2006; Reis, Fabian, & Miller 
[2010). While the source of black-body photons is strongly 
suspected to be the accretion disc, there is some controversy 
as to the origin and physical location of the power-law com- 
ponent, which could be produced by Compton scattering in 
a hot corona (Malzac & Belmont 2009), in a hot inner flow 
I ^Zdz iarski et al. 1998) resulting from the evaporation of the 
inner radii of the accretion disc (thus truncating the disc, see 
e.g. [Narayan fc Yi|1994| ), or at the base of the radio-emitting 
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jet that is observed during this state (Markoff, Nowak, &' 
[Wilms 2005). Notwithstanding the actual physical source of 
hard photons (i.e. E > 2 keV), there is general agreement 
that the hard photon-emitting region must be relatively cen- 
tral, within tens of gravitational radii from the central object 
in the bright hard state. 

Previously, most X-ray variability studies of the hard 
state were conducted using instruments with hard X-ray 
sensitivity, e.g. the Proportional Counter Array (PCA) on 
the Rossi X-ray Timing Explorer (RXTE), which limited 
the study of the variability to the hard, power-law compo- 
nent. Due to the anti-correlation of rms variability ampli- 
tude with spectral hardness as a source transitions to the 
disc-dominated soft state (e.g. Belloni et al. 2005), a pic- 
ture arose where the disc itself was intrinsically stable and 
variability is generated in the power-law emitting region, 
usually envisioned as an unstable hot inner flow (Churazov, 
Gilfanov , El Revnivtsevp601| |Done, GierSnski, &: Kubota| 

|2Q07| ). 

However, recent work (jWilkinson & Uttley|2009(|Uttley| 

'et al. 2011 ) has used the soft X-ray response and timing ca- 
pability of the EPIC-pn instrument on board XMM-Newton 
to carry out the first spectral-timing studies of the disc in 
the hard state. This work has shown strong evidence, in the 
form of variability ('covariance') spectra and time-lags, that 
the disc is intrinsically variable, and that disc variations in 
fact precede (and likely drive) the power-law variations, at 
least on time-scales of seconds and longer. The role of the 
disc as the driver of variability can be explained by invoking 
the presence of mass-accretion rate fluctuations that origi- 
nate in the disc at particular radii and propagate inwards. 
Such fluctuations vary at the local viscous time-scale corre- 
sponding to the radius where they are produced ( Lyubarskii 
1997; Kotov, Churazov, & Gilfanov 2001; Arevalo & Uttley 
2006). Because the m variations at inner radii are modulated 
by the outer rh via a multiplicative process in the context of 
this model, variability is seen at all frequencies that corre- 
spond to the viscous time-scales of the radii generating the 
propagating signals. Once the perturbations reach the inner 
regions of the system, these modulate the power-law emis- 
sion that accounts for most of the hard X-ray flux. A direct 
implication of this model is that the variations in power-law 
X-ray emission must be observed after the disc black-body 
variations with a time-delay scaling with the viscous travel- 
time between the radii where the fluctuations originate and 
the power-law emitting region. 

So far, evidence for disc-driven variability has been 
seen in XMM-Newton EPIC-pn observations of three hard 
state sources: GX 339-4, Cyg X-1 and SWIFT J1753.5- 
0127 (Uttley et al. 2011). These data were obtained during 
relatively long-lived hard states at moderate luminosities, 
around 1 per cent of the Eddington limit (with some un- 
certainty given the uncertain distances to GX 339-4 and 
SWIFT J1753. 5-0127 as well as some uncertainty on their 
masses). More luminous hard states, corresponding to out- 
burst rises or 'failed' state transitions, have previously not 
been studied with detailed spectral-timing into the soft X- 
ray band, but could contain important information about 
changes in the disc variability as the accretion rate increases 
and the disc emission strengthens. In this paper, we carry 
out the first detailed spectral-timing study of a luminous 
hard state of SWIFT J1753. 5-0127, a transient black hole 
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Figure 1. Swift/BAT lightcurve of SWIFT J1753. 5-0127, tak- 
ing weekly flux averages. The 2006 (#1) and 2009 (#2) XMM- 
Newton observations of the source are marked in red and orange, 
respectively. 

candidate which was discovered on 2005 May 30 (Palmer^ 
al. 2005). Its ~ 3.2 hr period derived from optical lightcurves 
( |Zurita et alp OOS Durant et al. 2009) makes it the BHXRB 
with the second-shortest period known to date (Kuulkers 
et al. 2012). Since the start of the outburst, it has under- 
gone a transition to and from the hard-intermediate state. It 
has never completed the outburst cycle toward softer states 
nor has it gone into quiescence. The long-term Swift/BAT 
lightcurve of SWIFT J1753. 5-0127 is shown in Figure 1 and 
shows the epochs of the two XMM-Newton EPIC-pn obser- 
vations. The first, obtained early in 2006, corresponds to 
a relatively faint hard state and shows clear evidence for 
disc-driven variability (Wilkinson & Uttley 2009; Uttley et| 
al. 2011). The second observation was triggered by us in 
response to a brightening of the hard state in September 
2009, and enables a comparison of a more luminous hard 
state with the lower-luminosity states studied to date. 

In Section 2 we describe our data reduction, and in Sec- 
tion 3 we carry out a detailed spectral-timing study of the 
2009 data and compare our results with the 2006 observa- 
tion. In Section 4 we summarise our key results and interpret 
them in terms of the disc fluctuation model which can ex- 
plain the data obtained at lower luminosities. 



2 DATA REDUCTION 

SWIFT J1753. 5-0127 was observed with XMM-Newton on 
2009 September 29 (ObsID 0605610301) and is shown as Obs 
#2 on Fig. 1, together with the 2006 observation (Obs #1, 
ObsID 0311590901) that we use later for comparison. We 
have extracted the data from the EPIC-pn camera in Timing 
Mode (Kuster et al. 2002) using the SAS 10.0.0 tool epproc 
and the most recent Current Calibration Files (CCFs). We 
have used evselect for flltering and extracted only single 
and double events (PATTERN <= 4) during intervals of low, 
quiet background. The column filters RAWX in [27:50] was 
chosen for the 2006 observation, and the selection RAWX in 
[27:35] II RAWX in [38:50] was adopted for the 2009, 
thus excising the central two columns to mitigate pile-up 
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effects. Although an excision of further columns results in a 
further decrease of pile-up, this effect does not change any of 
the results shown in this paper. Due to telemetry drop-outs, 
the total useful exposure are ~ 35.5 ks (2006) and ~ 13.3 ks 
(2009). 

The tool evselect was used to extract the spectrum 
and the response files were extracted with rmfgen and 
arfgen. Power spectra, cross-coherences, phase- and time- 
lags were extracted using data within the GTI intervals fol- 
lowing |Nowak et ahj (|l999| and |Vaughan fc Nowak| p997| | . 
We use a time bin of 0.005965 s and segment size of 16384 
and 4096 bins/segment (for the 2006 and 2009 observations, 
respectively) . 

We omit any segments with gaps between successive 
events longer than 0.1 s. These gaps are likely to be due to 
telemetry drop-outs that are not always included in the GTI 
files. Covariance (Wilkinson & Uttley 2009) and rms spec- 
tra (Revnivtsev, Gilfanov, & Churazov 1999) as well as the 
energy-dependent phase- and time-lags were also extracted 
for different frequency ranges as discussed in Section 3.3'| 

Given the Fourier transform of the signal in the i-th 
energy channel (or energy band) Siivi), the cross-spectrum 
between two bands 1 and 2 is defined as: 



{Cl,2(/^i)) = {Sl{vi)S2{u,)) 



(1) 



where the averages are performed over a number of inde- 
pendent segments and frequency bins m. Energy-dependent 
products such as phase- and time-lags as well as covari- 
ance spectra can be obtained from the above quantity by 
defining a reference band for each channel/band Sjiyi) = 
X^i'^^i Si' iyi)- In reality, the reference band can be fine-tuned 
to incorporate a given set of channels (e.g. a soft or a hard 
band), as we will show in the Results section. With this 
definition, the cross-spectrum for the i-th channel becomes: 
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Figure 2. 2009 spectrum (black) and best-fitting model (red 
solid line) using the model tbnew(l) * (diskbb(l)+iithconip(l) ) . 
The Ax^ and ratio residuals are shown in the lower panels. The 
disc black-body component is shown in blue (absorbed is shown 
dash-dotted, unabsorbed is shown dashed), whereas the Comp- 
tonisation component is depicted in green (absorbed is shown 
dash-dot-dotted, unabsorbed is shown dashed). The spectrum be- 
low 0.7 keV (left side of the vertical line) is shown but not used 
for the fit. 



{C.{vi)) = {Si{ui)S.{vi)) = I^S'Aj.i)S.{ui)^ (2) 

The energy-dependent phase-lag is then given by 
4)i{vi) = arg [{Ci(vi))] and its time-lag is Ti{vi) = 

(j)i{vi)/2-KUi. 

The covariance spectrum in the frequency range vi-^ to 
1^12 can be defined as ( Wilkinson!l|2011[|: 



CoY,{ui) = ^ , (3) 

where Avi is the width of the Z-th frequency bin and Bi{vi) 
denotes the bias caused by the contribution of Poisson power 
Ni , which is subtracted from the noiseless modulus-squared 
value of the cross-spectrum (see e.g. |Vaughan fc Nowakj 
,1997): 



BM) H{\SM)?){\S^{n)\') + {\SM)?){Nr{v:)) 
+ {N-,{,.i)){m{n)))/mi 



(4) 



The spectral analysis was performed in ISIS 1.6.2-1 
pfouck & Denicola 2000), using a spectrum where a min- 
imum of 3 channels per bin are grouped together, after 



adding 1% systematic error to account for instrumental un- 
certainties in the calibration of the instrument (Wilkinson] 
|&Uttley|2009|. 



3 SPECTRAL-TIMING ANALYSIS AND 
RESULTS 

3.1 Spectral analysis 

We first apply a simple phenomenological model to the en- 
ergy spectrum in order to understand qualitatively the con- 
tributions of the key continuum components, which will in- 
form our choice of energy-bands for spectral-timing analysis. 
A cursory examination of the unfolded spectrum reveals a 
clear soft X-ray excess above the usual power-law emission, 
consistent with a stronger disc black-body component than 
is seen at lower luminosities in the hard state. Therefore 
we consider a model consisting of a multi-colour disc black- 
body emission component plus a thermal-Comptonisation 
component (nthcomp, hereinafter the 'power-law' compo- 
nent given the similarity in shape at the hard energies), 
both absorbed by neutral Galactic absorption with the 
model tbnew ^ (Wilms et al.[,in prep.), assuming the abun- 
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Figure 3. Same as Fig.j2j, for the 2006 observation 



dances in Wilms, Allen, & McCray (2000) and the cross- 
sections in Verner et al. (1996). In ISIS, the model cor- 
responds to tbnew(l) * (diskbb(l)+nthcomp(l) ) . In addi- 
tion, we fix the neutral hydrogen column density A^'^r value 
to 0.194 X 10^^ cm"^ (Wilkinson & Uttley 2009) and the 
thermal-Comptonisation electron temperature kTe to 53 
keV (Reynolds et al. 2010). The seed photon temperature for 
Comptonisation is tied to the disc black-body inner disc tem- 
perature. The model is fitted only to the 0.7-10 keV energy 
range, following the recommendations of the EPIC Consor- 
tium not to consider the data below 0.6 keV (Guainazzi 
^llp , however we plot the entire useful energy range cov- 
ered by the data for completeness and to show clearly the 
predicted disc contribution. 

The best-fitting parameters for the model are shown 
in the first column of Table 1 and the spectrum together 
with the best-fitting model is shown in Figure 2 Although 
the fitting procedure yields a high x^/d.o.f. ~ 2.6 (possibly 
because of underestimated instrumental systematic errors), 
our aim is not to provide an accurate spectral fit but rather 
to understand what emission components are the fundamen- 
tal spectral sources of photons. Therefore, we do not include 
the non-primary spectral components in the fit, such as the 
disc reflection component. 

We use the same model to re-fit the EPIC-pn spec- 
trum of the 2006 hard-state observation presented in Miller, 
Homan, fc Miniutti (2006) and Wilkinson & Uttley (2009), 



which we will later use for a comparison with the 2009 data. 
The fitted parameters are shown in Table 1 and the model 
fit is plotted in Fig. 3 In this case we find a harder power- 
law (r = 1.620 ±0.005) and a much weaker disc black-body 
component with a lower temperature (feTin = 0.212^Q Q^g) 
and normalisation only ~ 3 percent of the value from 2009. 



Table 1. Best-fitting values for the model 
tbnew(l) * (diskbb(l)+nthcomp(l) ) fitted over the range 
0.7 — 10.0 keV. The disc black-body normalisation N^i^y^ is in 
units of [(ri,i/km)/(Z)/10kpc)]^ cos 9 (where rjn is the 'apparent' 
disc inner radius that is also dependent on the inclination 9 of 
the disc, see Kubota et al. (1998)), whereas the Comptonisation 
normalisation is in units of photons cm~^ s~^ keV~^ at 1 keV. 
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Value (2009) Value (2006) 



Nh [10^2 cm-2] 


0.194 (fixed) 


0.194 (fixed) 


fcTin [keV] 


0.248 ±0.001 


r, 9+0.020 


fcTc [keV] 


53 (fixed) 


53 (fixed) 


A^dbb 


314901^^;; 


1050t^™ 


F 
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1.620 ±0.005 
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Figure 4. Power spectral densities (PSDs) for the 2009 obser- 
vation in the energy bands 0.3 — 0.6 keV (black triangles), 0.6 - 
1.0 keV (red squares), 1.0 - 2.5 keV (blue diamonds), 2.5 - 4.0 
keV (green plus signs) and 4.0 — 10.0 keV (magenta times signs) 
showing two broad peaks whose relative amplitude varies as a 
function of the energy band chosen. 



3.2 2009 data: frequency-dependent 
spectral-timing properties 

The presence of both a disc black-body component and a 
Comptonisation component that dominate different energy 
bands covered by EPIC-pn, as seen in Figure 2, suggests 
that the variability properties of the source may differ where 
the photon flux is dominated by either the disc black-body 
or the power-law components. Hence it is convenient to ex- 
tract light curves for the energy ranges 0.3 - 0.6 keV, 0.6 - 
1.0 keV, 2.5 - 4.0 keV and 4.0 - 10.0 keV, hereinafter ul- 
trasoft, soft, intermediate and hard bands, respectively. For 
completeness we also extracted data from the the 1.0 - 2.5 
keV band, although since this band shows no unique prop- 
erties compared to the other bands, we do not show it in our 
later spectral-timing results, but we do include it in show- 
ing the energy-dependent power-spectral densities (PSDs), 
which are shown in Fig. 4, As shown in Fig. 2, the ultrasoft 
and soft bands are dominated by the disc black-body com- 
ponent, whereas the intermediate and hard bands cover the 
power-law component. The range 1.0 - 2.5 keV represents 
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Figure 5. Coherence, time- and phase-lags for the 2009 observation, using different combinations of the bands 0.3 — 0.6 keV, 0.6 - I.O 
keV and 4.0 - 10.0 keV. A positive lag indicates that the emission in the harder band lags the emission in the softer band. The shaded 
regions indicate the frequency intervals used to average the energy-dependent products studied in Section 3.3 



the variability power for the photons emitted by a combina- 
tion of the black-body and power-law components. 

In addition, the ultrasoft, soft, intermediate and hard 
bands are also used to extract the coherence, phase- and 
time-lags, using the combinations ultrasoft-soft, ultrasoft- 
intermediate, ultrasoft-hard, soft-intermediate, soft-hard 
and intermediate-hard, shown in Fig. 5, The properties of 
coherence and lags with respect to the 1.0 - 2.5 keV band 
show an intermediate behaviour compared to other bands 
and are not shown here. 

3.2.1 Amplitude of variabihty 

The power spectral densities in Fig. 4 show evidence for 
two broad noise components, peaking at ~ 0.4 and ~ 2 
Hz (hereinafter, the 'low-frequency' and 'high-frequency' 
peaks respectively). The relative amplitude between the two 
peaks changes, depending on the energy band considered. 
At the soft and ultrasoft energies where the disc black-body 
emission dominates, the low- frequency peak shows a larger 



peak-amplitude of variability than the high-frequency peak. 
Conversely, for the power-law-dominated intermediate and 
hard energies the high-frequency peak-amplitude becomes 
stronger with respect to the low-frequency peak and is con- 
sistent between the two bands, while the 1.0 - 2.5 keV band 
shows an intermediate behaviour between the softer and the 
higher-energy bands. Besides the energy-dependence of the 
PSD shape, there is also a clear change in normalisation, 
with the ultrasoft band showing the weakest variability at all 
frequencies, followed by the soft band and then the higher- 
energy bands. This behaviour contrasts with that seen in the 
2006 data by Wilkinson & Uttley (2009), where the softer 
0.5-1 keV band showed larger- amplitude variability than the 
2-10 keV band at low frequencies. 

3.2.2 Coherence and frequency- dependent spectra 

The coherence is an indicator of the degree of linear cor- 
relation between variations in the emission in two different 
bands (Vaughan & Nowak 1997) at a particular frequency. 
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Full coherence (when the coherence equals unity at the given 
frequency) is reached when the variations in emission in one 
band can be reconstructed linearly from the variations in 
emission in the other band. The upper panels in Fig.'S'show 
the coherence computed for all possible combinations of the 
above bands. For all plots in the figure (including time- and 
phase-lags), the left panels show the frequency-dependent 
products computed between the pair of bands that cover 
the disc black-body region or the pair of bands which cover 
the power-law region, i.e. the disc vs. disc behaviour and 
the power-law vs. power-law behaviour. The remaining pan- 
els show the same products made from intercombinations of 
the disc and power-law bands, combining the intermediate 
and hard bands (middle and right panels, respectively) to- 
gether with the ultrasoft and soft bands that are dominated 
by the disc black-body emission. 

In a simple case where all of the emission in one spectral 
component varies linearly together at all energies, the coher- 
ence between the observed emission in two bands where that 
component dominates should approach unity. This is similar 
to the case observed with the coherence between the inter- 
mediate and hard bands that sample the power-law (C.i), 
as well as between the intermediate (C.ll) and hard bands 
(Cm) with respect to the soft band, up to ~ 0.6 Hz. This 
frequency roughly corresponds to the point where the low 
and high-frequency peaks in the PSD intersect (Fig. |4|). The 
coherence deviates slightly from unity for frequencies below 
~ 0.1 Hz between the hard and soft bands. 

A very different behaviour is observed between the ul- 
trasoft and the soft, intermediate and hard bands. In this 
case, unity coherence is reached at ~ 0.6 Hz, and decreases 
towards both lower and higher frequencies, reaching values of 
~ 0.5 in some cases. Interestingly, the ultrasoft variations at 
low frequencies are just as incoherent with the soft variations 
as they are with the hard or intermediate variations. Thus, 
at low frequencies, while the variations in the two power-law 
bands are coherent with each other, a significant fraction of 
the variations in the two disc bands are incoherent with each 
other. 



3.2.3 Time and phase lags 

Previous observations of time-lags in the hard state show 
that they roughly follow a frequency-dependence t{i') oc 



(Miyamoto & Kitamoto 1989 



[Miyamoto et al. 1992; 



[Crary et al.|1998||Nowak et al.|1999| |, albeit with detectable 
'steps' in the lag vs. frequency dependence where the data 

are good enough to discern them (Nowak et al. 1999). For the 
observation that we are studying, the frequency dependence 
does clearly vary as a function of energy bands chosen. For 
lags between the two disc-dominated bands (ultrasoft and 
soft) or between the intermediate and hard power-law bands, 
a relatively small amplitude is seen with a weak and smooth 
evolution of phase lag with Fourier-frequency. However, the 
lags between power-law and disc bands show a much more 
pronounced change with a significant 'step' (most clearly 
visible in the phase lags) above 1 Hz. Thus the strongest 
evolution of lags with frequency occurs when comparing the 
lags between the power-law and disc bands, with the power- 
law lagging the disc variations, as seen in other observations 
of hard state BHXRBs ) |Uttley et al.|2011 l. 



3.3 2009 data: energy-dependent lags and 
covariance spectra 

Given the complex behaviour in the coherence and lags 
shown in Fig. [5] we next pursue an energy-resolved analysis 
that can more clearly relate this behaviour to the different 
spectral components that contribute to the variability of the 

source. 

Wc extract cross-spectra in the frequency ranges 0.04 - 
0.20 Hz, 0.2 - 0.6 Hz, 0.6 - 1.0 Hz and 1.0 - 2.0 Hz (that 
correspond to each of the shaded regions in Fig. |5l that cor- 
respond to frequency ranges with comparable behaviour in 
terms of lags and coherence. The cross-spectra are computed 
between a broad, high signal-to-noise reference energy band, 
and each of 30 rebinned energy channels. Following |Wilkii>l 
Ison & Uttleyl (120091) and "Uttley et al.''(2011), the reference 
band light curve is determined separately for each channel 
since it must have the signal in the channel of interest sub- 
tracted in order to avoid correlated Poisson noise effects. In 
these previous works, the data were largely coherent at all 
frequencies, so the results were not sensitive to the choice 
of reference band. In the observation presented here, our 
spectral-timing products could vary with the choice of ref- 
erence band, because the lags and covariance spectra are 
effectively weighted by the coherence of each channel with 
the reference band. Therefore, to account for any effects due 
to the choice of reference band, we make cross-spectra for 
two re ference bands: 0.5 - 2.5 keV and 3.1 - 10.0 keV. Fol- 
lowing lUttley et al.| (|2011| , we make energy-dependent lag 
and covariance spectra for each of these two bands, plotted 
as the left and right columns in Fig. |6J 



3.3.1 Lag- energy spectra 

The time and phase-lag versus energy are shown in the top 
two sets of panels in Fig. 6 Because the lag difference be- 
tween two separate energy channels is relative to the choice 
of reference band, the offset on the j/-axis is not meaning- 
ful in this case. The main consideration is the relative lag 
between channels, which is plotted here so that more posi- 
tive values of lag are lagging smallcr/more-ncgativc values. 
We do not plot the time lags for the lowest-frequency range, 
since their large values would make it difficult to read the 
data for the other frequency ranges (these data can be seen 
in the phase lag plots, however). 

The dependence of the lags on energy can be approxi- 
mated as a 'broken' log-linear law, with a steeper slope be- 
low ~ 1.0 keV and a flatter dependence above that energy. 
This turning point is close to the energy where the photon 
flux from the Comptonisation component overcomes the disc 
black-body photon flux (Fig. 2 ) . In the frequency range 0.2 - 
0.6 Hz, the lag between the power-law and disc components 
reaches values up to ~ 0.2 s, reaching values close to zero as 
the frequency increases. This behaviour was also shown in 
Uttley et al. (2011) and is consistent with the disc variations 
leading (and probably driving) the power-law variations on 
time-scales of seconds. At frequencies of 1 - 2 Hz, the low- 
energy down-turn in lags is replaced by a small up-turn, 
also similar to what was observed for other observations of 
BHXRBs at similar frequencies ( jUttley et al.||2011| ). How- 
ever, it is interesting to note that in this frequency range, the 
log-linear shape appears also to break downwards at around 
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Figure 6. Time- and phase-lag spectra, covariance spectrum divided by average spectrum and covariance spectrum divided by covariance 
spectrum computed in the range 0.2 — 0.6 Hz for 2009 observation. The frequency ranges used are 0.04 - 0.20 Hz (black triangles), 0.2 
— 0.6 Hz (red squares), 0.6 - 1.0 (blue diamonds) and 1.0 — 2.0 Hz (green plus signs). The left and right plots show the products above 
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2 keV (this behaviour is seen most clearly in the phase lag 
plot), suggesting that there is a composite shape, perhaps 
consisting of leading higher-energy disc emission together 
with lagging lower-energy disc emission. It should be noted 
that the general behaviour in the lags is replicated using 
both choices of reference band, albeit with small differences 
between them. 

3.3.2 Covariance spectra 

In addition to the causal information shown above, we ex- 
tract covariance spectra from the cross-spectra (Wilkinson 
& Uttley 2009, Uttley et al. 2011) in order to quantify, as a 
function of energy, the amplitude of variable emission that 
is linearly-correlated with the reference band of choice (see 
Section 2). The covariance spectrum depends on the coher- 
ence between the reference band of choice and each energy 
channel. In the limit where the coherence is unity, the covari- 
ance spectrum is equal to the rms spectrum (e.g. "Revnivt-" 
sev, Gilfanov, & Churazov 1999), although with consider- 
ably smaller error-bars due to the effect of the broad refer- 
ence band acting as a matched filter on the flux variations 
in individual channels. 

The panels in the third row of Fig. 6 show the corre- 
sponding covariance spectra divided by the mean spectrum, 
which enables us to estimate the contribution of each com- 
ponent to the overall variability, analogous to the fractional 
rms spectrum. For energies above ~ 2 keV, the fractional 
covariance spectra show a roughly power-law slope that be- 
comes harder for increasing frequencies. At the lowest fre- 
quencies, the power-law slope is softer than the mean, while 
it reaches the slope of the average spectrum at the highest 
frequencies, 1.0 - 2.0 Hz. This effect is mild if only the first 
three frequency ranges are considered, and becomes more 
obvious for the 1.0 - 2.0 Hz range. 

Such an effect can also be seen in the lowest panels 
of Fig. 6, by computing the ratio between the covariance 
in the frequency bands 0.04 - 0.20 Hz, 0.6 - 1.0 Hz and 
1.0 - 2.0 Hz and the covariance in the range 0.2 - 0.6 Hz. 
In this case, the variation in power-law slope is plotted in- 
dependently of any constant component, which would not 
contribute to the covariance spectra used to normalise the 
data. Thus the effect cannot be due to any simple change 
in the ratio of a variable to constant spectral component 
which have different spectral slopes. Instead, the hardening 
of the power-law variations at higher frequencies is intrin- 
sic to the variable power-law itself. Similar behaviour has 
been seen in frequency-dependent rms spectra of hard state 
BHXRBs obtained from RXTE data (e.g. ||Revnivtsev, Gil-| 
ifanov, & Churazov 1999). 

At the disc-dominated energies (i.e. below ~ 1 keV), 
the fractional covariance is up to a factor ~ 5 weaker than 
at the power-law-dominated energies, suggesting that, al- 
though the disc variability is apparently driving the power- 
law variability (as seen in the causality argument provided 
with the lags above, and in Uttley et al. 2011), the disc 
emission is less variable than the power-law emission. 

In Fig. 5 we showed that the coherence clearly drops 
below unity at low frequencies in any of the combinations of 
energy bands studied where the ultrasoft band was consid- 
ered. As in the previous cases shown above, this drop can 
be quantified in an energy-dependent manner by comparing 
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Figure 8. Power spectral densities as in Fig. 4, for the 2006 
observation 

the rms and covariance spectra. The fractional covariance 
spectra shown in Fig. 6 appear to flatten off at the low- 
est energies in the intermediate frequency bands, while in 
the lowest-frequency band (where coherence is low) it con- 
tinues to drop. We can show that this effect is due to the 
low coherence of variations of the softest photons by com- 
paring the fractional covariance spectra with the fractional 
rms spectra, shown in Fig. 71 Where the fractional covari- 
ance and rms are the same, the coherence is consistent with 
unity. However, at low frequencies, the fractional covariance 
is clearly lower than the fractional rms below ~ 0.5 keV, 
indicating sub-unity coherence at these energies, consistent 
with our findings from the frequency-dependent coherence. 

On the other hand, above 0.5 keV the fractional rms 
and covariance are similar for the three lowest frequency 
ranges, indicating that the higher-energy part of the disc 
emission varies coherently with the power-law emission on 
time-scales longer than a second. However, as also expected 
from the frequency-dependent coherence measurements, the 
coherence of this soft-band emission with the higher energy 
emission also drops in the 1-2 Hz frequency bin. This drop 
manifests itself as a drop in fractional covariance compared 
to fractional rms, which is seen at low energies for the 3.1- 
10 keV reference band, and at high energies for the 0.5- 
2.0 keV reference band (since this band contains the disc 
emission). 

3.4 Comparison with the spectral-timing 
properties in the 2006 hard state 

The 2006 observation was obtained during a hard state that 
was a factor of ~ 3.4 fainter than the 2009 observation, as 
seen in the Swift/BAT lightcurve (red arrow in Fig. 1). The 
luminosity difference is even greater at the lower energies 
covered by XMM-Newton, since the 2006 data show a much 
harder power-law component (F = 1.620±0.005) which dom- 
inates the SED, since the disc component is much weaker 
(see Figs. 2 and 3 and Table 1 for a comparison). Due to 
the relatively low signal-to-noise in the 2006 data compared 
to 2009, we limit our comparison of the data to a few key 
aspects (such as the PSD and coherence) and the two lower 
frequency ranges for the energy-dependent comparison of 
lags and covariance. 
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Figure 7. Fractional covariance (black triangles) and fractional rms (red squares) in the 2009 observation for the frequency ranges 0.04 

- 0.20 Hz, 0.2 - 0.6 Hz, 0.6 - 1.0 Hz and 1.0 - 2.0 Hz. The left and right plots show the products above using the 0.5 - 2.5 keV and 3.1 

— 10.0 keV reference bands, respectively. 




0.1 

Frequency [Hz 



Figure 9. Coherence values for the 2006 observation, using the same energy bands as in Fig. [s] 



The power spectral densities for the 2006 observation 
are shown in Fig. 8 In this observation, the data can be ap- 
proximated by a power-law down to ~ 0.04 Hz with a similar 
slope and amplitude for all energy ranges considered, as op- 
posed to the much clearer energy dependence of the power 



amplitude of the two peaks observed in 2009, whose relative 
amplitudes also changed with energy. The lower-frequency 
peak in the PSD seen in 2006 is consistent with the correla- 
tion of PSD changes with spectral evolution seen in moder- 
ate to high luminosity hard states of other BHXRBs, with 
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lower-frequency features seen when the spectrum is harder 
(e.g. Belloni et al. 2005, also seen in Bock et al. 2011). 

In the 2009 observation, an important property of the 
coherence at frequencies below ~ 0.6 Hz is a drop down to 
coherence values ~ 0.5 for each combination of bands that 
involved the ultrasoft band. The corresponding coherence 
plot for 2006 (Fig. 9) shows a value much closer to unity 
down to the lowest-frequencies probed. 

Fig. 10 shows the energy-dependent lags and frac- 
tional covariance spectra for the 2006 data for the first two 
frequency-ranges (0.04 - 0.20 Hz and 0.2 - 0.6 Hz), using 
the reference band 0.5 - 2.5 keV. The harder reference band 
gives consistent but noisier results (as expected, since the 
coherence is high at these frequencies). These plots show 
behaviour consistent with that previously reported for this 
observation by Wilkinson & Uttley (2009) and Uttley et 
|al. (2011). The lowest-frequency lags show the characteris- 
tic drop below 1 keV that can be associated with the disc 
driving the variability at harder energies, consistent with 
the behaviour also seen in 2009 (and implying that there 
is intrinsic disc variability in both observations). However, 
the covariance spectra show a 'soft excess' at the lowest fre- 
quencies, implying that there is extra disc variability rela- 
tive to the power-law variations at these frequencies. This 
behaviour contrasts with that seen in 2009 where the disc 
variations at all frequencies appear to be suppressed when 
compared with the power-law variations. 

3.5 The emissivity profiles 

The spectral fits to the 2009 data show much stronger disc 
emission than is seen in 2006 (Table 1). Although the in- 
ferred disc temperatures are relatively similar, the normal- 
isations (which are proportional to the square of the inner 
disc radius) are very different. 

Taken at face value, the disc black-body normalisation 
in 2009 implies an inner disc radius 5-6 times larger than 
that in 2006, which is at odds with the standard interpre- 
tation of the other phenomenological changes between the 
two data sets, which assumes that the disc becomes less trun- 
cated as PSD characteristic frequencies increase and the disc 
emission becomes stronger. However, to interpret changes in 
disc black-body normalisation in terms of changes in inner 
disc radius, we must assume in that both data sets the disc 
emissivity profile is the same. Since the X-ray luminosity 
in 2006 is dominated by the power-law emission, it is likely 
that much of the disc emission is produced by X-ray heating 
of the disc by the power-law, which could produce a much 
more centrally peaked (and apparently 'smaller') disc black- 
body component (e.g. see Gierlihski, Done, & Page 2008). In 
2009, the disc emission is comparable to or even exceeds that 
from the power-law and most of the emission is likely to be 
intrinsic to the disc, leading to a more 'standard' emissivity 
profile. Since the peak energy of the disc black-body emis- 
sion is comparable to the low energy end of the spectrum 
included in the fit, distinguishing among different emissiv- 
ity profiles from the spectral fits would be an arduous task. 
Systematic errors deriving from spectral calibration do not, 
however, impact on the variability analysis performed here. 

Assuming that the disc emissivity profile in 2009 has 
the standard form (in terms of bolometric emission) 

and hence the disc temperature scales with radius as J?""^^*, 
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Figure 10. Time- and phase-lag spectra and covariance spectrum 
divided by average spectrum for 2006 observation, using the 0.5 
- 2.5 keV reference band. The frequency ranges used are 0.04 - 
0.20 Hz (black) and 0.2 - 0.6 Hz (red). 



we can consider the dependence of cumulative disc emission 
on radius.^ 



^ In their optical/UV-emitting regions, the temperature profiles 
of irradiated discs should be flatter than those of disks which are 
purely viscously heated, if the additional heating increases the 
scale-height of the disc (e.g. Vrtilek et al. 1990) or if the disc 
is signiflcantly warped. However, the inner, soft X-ray emitting 
regions of discs considered in this paper are likely to be heated 
primarily by viscous heating, with, in the case of the 2006 data, 
an additional contribution from irradiation which is concentrated 
in the centre (i.e. assuming a compact corona). Outside this in- 
nermost region of intense irradiation, the discs are likely to have 
temperature proflles closer to those predicted from the viscous 



For comparison and allowing for the likely possibility 
of substantial disc-heating by the power-law in 2006, we fit 
the 2006 spectrum with a combined single black-body (rep- 
resenting a hot inner ring heated by the central power-law 
emission) and a standard disc black-body component. We 
force the disc black-body to have the same normalisation 
(i.e. the same assumed inner radius) as the disc in 2009 
(note that due to degeneracy in the model the inner radius 
in 2006 may be even larger, as expected if the disc is more 
truncated). For this composite model the disc black-body 
temperature drops to 0.11 keV while the inner single black- 
body temperature is high at 0.31 keV, but corresponds to 
only a small emitting area at the inner edge of the disc. The 
best-fitting model parameters give x^/d.o.f.= 565/623, pro- 
viding an even better fit to the data than the single disc 
black-body model. The resulting cumulative emission pro- 
files are shown in Fig. 11 for the soft and ultrasoft energy 
bands (after allowing for the efi'ects of the instrumental re- 
sponse on each band). 



4 DISCUSSION 

To summarise, our main observational findings are: 

1. The 2009 XMM-Newton observation of 
SWIFT J1753. 5-0127 took place when the source was 
in a significantly more luminous hard state (~ 6 times more 
luminous in the 0.5-10.0 keV band) than in 2006. Corre- 
spondingly, the power-law spectral component is softer than 
in 2006, and importantly the soft excess associated with 
the disc black-body is significantly stronger. The inferred 
black-body temperature in 2009 is only marginally higher 
than in 2006, but the 2009 disc black-body normalisation is 
a factor ~ 30 larger than in 2006. The larger inferred inner 
disc radius in 2009 conflicts with standard expectations for 
reduced disc truncation as the disc luminosity increases. 
However, the data can be reconciled with a non-changing 
or even a smaller inner radius in 2009, if the intrinsic disc 
temperature in 2006 is significantly lower than inferred 
from the single disc blackbody fit and there is additional 
blackbody emission from a hot inner ring, due to X-ray 
heating by the relatively stronger power-law in 2006. 

2. The PSDs in 2009 are band-limited, and can be inter- 
preted in terms of two broad frequency components which 
peak around 0.4 Hz and 2 Hz respectively. 

3. The 2009 PSD shape and amplitude evolves with en- 
ergy. The lowest-amplitude variability over all frequencies is 
seen in the softest 0.3-0.6 keV band, followed by the 0.6- 
1.0 keV band. The amplitude of the lowest-frequency com- 
ponent is largest in the 1.0-2.5 keV band before dropping 
slightly at higher energies, while the amplitude of the high 
frequency component increases with energy. 

4. The covariance and rms spectra suggest that the 
energy-dependent PSD behaviour is linked to two effects. 
Firstly, there is a fall-off in fractional rms (seen in all fre- 
quency ranges) at low-energies where the disc black-body 
dominates the emission. Secondly, the variable power-law 



heating models, only flattening on larger scales where soft X-rays 
are not emitted. 
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Figure 11. Normalised cumulative emission function in the ul- 
trasoft (solid black) and soft bands (dashed red) as a function of 
disc radius. The top panel corresponds to the 2006 data, assuming 
a composite disc black-body and inner hotter black-body associ- 
ated with X-ray heating by the power-law emission (see text for 
details). The bottom panel corresponds to the 2009 data, assum- 
ing a dominant intrinsic disc black-body profile. 



emission component becomes harder at higher temporal fre- 
quencies. Neither of these effects are seen in the 2006 data. 
In fact, in 2006 the opposite is seen at low energies, where 
the fractional rms increases in the energy range where the 
disc contributes to the spectrum. 

5. Frequency-dependent lag measurements show that in 
general, the harder photons lag softer ones. Within the en- 
ergy bands dominated by the power-law (2.5-4 keV and 4- 
10 keV) and the disc (0.3-0.6 keV and 0.6-1 keV) the lags are 
relatively short and the phase-lags evolve only smoothly and 
weakly (if at all) with frequency. However, the lags between 
power-law and disc bands are larger and show significant 
structure in the form of 'steps', with a significant increase 
in phase lag with frequency. 

6. Lag-energy spectra show similar behaviour to that re- 
ported by us (Uttley et al. 2011) in other hard state observa- 
tions of BHXRBs (including SWIFT J1753. 5-0127 in 2006). 
The sharpest change in the lags is seen at low energies where 
the disc dominates the spectrum. The low-energy shape of 
the lag spectrum depends on the frequency probed: the high- 
est frequencies show evidence for a shift in the break in the 
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lag spectrum to higher energies, together with an upturn in 
the lowest energy bins (similar to that seen in GX 339-4, 
[Uttley et al.|2011| |. 

7. The coherence between the bands dominated by the 
power-law remains high over a broad frequency range. How- 
ever, between the power-law and disc-dominated bands, and 
between the two disc-dominated bands, coherence drops sig- 
nificantly at frequencies below 0.5 Hz (where the softest disc 
band shows the least coherence with the other bands) as well 
as above 1 Hz (where both disc bands show low coherence 
with respect to the power-law bands). Comparison of the 
covariance and rms spectra confirm that low-frequency vari- 
ations of the softest photons, below 0.5 kcV, arc only weakly 
correlated with variations at higher energies. In contrast, the 
2006 data are relatively coherent at low frequencies. 

The spectral-timing behaviour of SWIFT J1753.5 0127 in 
2009 is complex but we will attempt to interpret it here 
in terms of the physical picture of mass-accretion rate fluc- 
tuations arising in and propagating through the accretion 
disc, suggested by Wilkinson & Uttley, (2009) and^ttley et 
al. 1 2011[ ) in order to explain the covariance and time-lag 
behaviour of hard state BHXRBs. Since we see similar lag 
behaviour in the SWIFT J1753. 5-0127 2009 data, it is use- 
ful to understand whether this model can be extended to 
explore the more complex spectral-timing behaviour seen in 
this datasct. In doing so, we will focus primarily on explain- 
ing the differences between the spectral-timing properties 
of the 2009 data we present here and those from 2006 and 
observations of other BHXRBs reported in Wilkinson & Ut- 
|tley_ (2009 ) and Uttley et al. (2011). Comprehensive testing 
of the propagating fluctuations model against these data is 
beyond the scope of this work, so here we simply sketch 
how the key observational results might be explained by the 
model. 



4.1 Implications from the propagating 
fluctuations model 

The key feature of the propagating fluctuations model used 
to explain the soft X-ray behaviour of other hard state 
BHXRBs is that the disc plays an important role in gener- 
ating and carrying the fluctuations (at least on fluctuation 
time-scales > 1 s). In Subsection 3.5 we assumed that the 



intrinsic disc emissivity profile is the same in both observa- 
tions (oc R~^^'^). Under this scenario, we showed that the 
disc emission in 2006 may be explained if the overall disc 
emission in 2006 is the result of the extended emission with 
a 'standard' emissivity profile, plus a thermal emission com- 
ponent that is caused by significant heating of the disc by 
the illuminating power-law photons. 

This additional inner-disc emission produces up to ~ 75 
per cent of the disc emission in the soft band (Fig.[Tll). Due 
to the lower disc temperature and the relatively hot inner 
ring, the disc emission in 2006 is much more centrally con- 
centrated than in 2009. Furthermore, the ultrasoft emission 
is significantly more extended across the disc than the soft 
emission. 

We can use these differences to explain some of the in- 
teresting spectral-timing behaviour in the data. Note that 
the arguments that follow apply independently of the actual 
inner disc radius rin, since the emission profile in each band 



is set by the disc temperature at rin (which is set by the 
spectral fit results) and the relative change in radius. 

First, let us consider the 2009 data and suppose that 
a fluctuation in the disc arises at a relative small radius 
R ~ 1.5 Tin, and modulates the disc emission as it propa- 
gates inwards. Since only ~40 per cent of the ultrasoft flux 
is contained within this radius, the observed variability am- 
plitude will be suppressed by a factor ~ 2.5 compared to the 
amplitude of the intrinsic fluctuation. In contrast, the soft 
band fluctuation will be suppressed by only a factor ~ 1.5, 
while the power-law emission, if it originates within the disc 
truncation radius, will not be suppressed at all. This sim- 
ple i^icturc can explain the energy dependence of the PSD 
amplitude in 2009 and the corresponding drop in rms and 
covariance spectra at low energies. In this scenario, the same 
behaviour may not seen in other hard-state XMM-Newton 
observations of BHXRBs or in the 2006 data because of 
two effects. Firstly, since the PSDs show lower characteristic 
frequencies in these other observations, the fluctuations in 
accretion rate may arise at a larger radius, so that a sigiufi- 
cantly larger fraction of the disc emission can be modulated. 
Secondly, as can be seen in the 2006 emission profile, the 
combination of X-ray heating of an inner ring with a cooler 
disc at larger radii will lead to more concentrated emission 
in both disc bands, so that all or most of the flux will be 
modulated by a fluctuation arising at i? ~ 1.5 nn or larger. 
Thus we do not see suppression of the PSD amplitude at 
soft X-ray energies. 

4.1.1 Coherence and frequency-dependent spectra 

We next consider the drop in coherence at low frequencies. 
The natural explanation for this behaviour is that there is 
some component of accretion fluctuations in the disc which 
are unable to propagate very far inwards and therefore ap- 
pear as variations in the extended disc emission (in the ul- 
trasoft band and to a lesser extent the soft band) but do 
not produce correlated variations in the power-law emis- 
sion. One possibility is that these variations are generated at 
larger radii than the spectrally-coherent fluctuations, so that 
they are viscously damped before they are able to propagate 
to the innermost part of the disc. Again, we may see this 
effect in the 2009 data because of the strong intrinsic disc 
emission and extended emissivity profile. Disc emission vari- 
ations driven by power-law heating of the disc will be much 
more coherent with the power-law variations. Furthermore, 
the lower-frequencies seen in the PSD of the 2006 data could 
imply that signals can effectively propagate from larger radii 
in the disc and thus modulate a significant fraction of the 
intrinsic disc emission coherently. 

The higher-frequency drop in coherence corresponds 
roughly to the overlap between the low and high-frequency 
PSD components, and could be produced by a similar effect 
to that described above. For example, if the higher-frequency 
signal corresponds to a second, even smaller radius in the 
disc which generates propagating accretion fluctuations, the 
signals generated at the same frequencies at larger radii may 
not propagate effectively through the disc, again producing 
fluctuations that are not seen in the power-law component. 
It is important to note however that the relatively small radii 
which we suggest as the origin of fluctuations (e.g. within 
1.5 Tin) also imply that the power-law emitting region must 
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be compact. Otherwise, seed photons from the varying re- 
gions in the disc would upscatter to produce variations of 
the power-law that are coherent with the disc variations. 
Similarly, we expect some component of disc emission to be 
produced by X-ray heating by the power-law and hence be 
coherent with the power-law variations. These components 
are likely to exist at some level in the variable X-ray signal, 
but they must be relatively weak, at least on the time-scales 
where the coherence is low. 

Some level of interaction between the disc jjlioton vari- 
ations and the central power-law emitting region could po- 
tentially explain the observed hardening of the covariance 
and rms spectra at higher frequencies (and correspondingly, 
the peak in low-frequency power in the 1-2.5 kcV energy 
range). The amplitude of soft-photon variations is largest 
in the low-frequency peak of the PSD. Compton cooling of 
the power-law emitting region by these soft photons could 
lead to a significant pivoting of the power-law in response 
to these variations, such that at low frequencies the frac- 
tional amplitude of power-law variations increases towards 
lower energies. This effect then has the converse implication 
that at high frequencies the variable power-law appears to 
harden. 



around 1 Hz). It is important to note however that the rise 
in the time lag seems to coincide with the frequency range 
where the coherence peaks around unity. Incoherent signals 
which are caused by a mixture of uncorrelated variable com- 
ponents will show lags which are a mixture of the lags of each 
separate component. Thus, the rise in lag could be explained 
if the component causing the low-frequency drop in coher- 
ence contributes a small lag, reducing the total observed 
lag, but at higher frequencies of ~ 0.5 Hz this component 
disappears from the variable signal and the lag suddenly in- 
creases as a result. One possibility is that the incoherent 
fluctuations in the disc also modulate a small disc-corona 
component which is physically separate (and also spectrally 
distinct from) the inner power-law emitting region, which 
may be linked to a hot inner flow or the base of a jet. The 
interband lags of such a component would then be relatively 
small, and would give way to the larger lags (due to propa- 
gation effects) once the coherent signal becomes dominant. 
If the lags of the coherent signal themselves show a depen- 
dence on frequency, (e.g. a weak increase in phase lag with 
frequency as observed in other hard state BHXRBs) a com- 
plex pattern of rises and falls could be produced. 



4.1.2 Lags 

Finally, we consider the lag behaviour. The relatively small 
lags and smooth dependence on frctjucncy seen between the 
two power-law bands and also the two disc bands can be sim- 
ply explained because here we are seeing the same spectral 
components. Fluctuations propagating through the disc will 
primarily modulate the ultrasoft emission first, but there is 
significant overlap between soft and ultrasoft emission across 
the disc, so we expect lags to be relatively small. The same is 
true for the power-law component - lags should be small be- 
tween energies where the power-law component dominates. 
The largest lags are then expected between the disc and the 
power-law components, assuming that these are physically 
separated, i.e. the power-law emission primarily originates 
within the disc truncation radius. 

We must still explain the clear structure in the lag vs. 
frequency dependence of the power-law relative to the disc- 
dominated bands, i.e. the apparent 'steps'. Similar steps can 
be seen in other hard state data from RXTE, but between 
power-law dominated energies (Nowak et al. 1999). In the 
SWIFT J1753. 5-0127 2009 data, the effect is particularly 
strong: the time lags actually increase with frequency above 
~ 0.5 Hz before dropping above ~ 1.5 Hz. |Nowak| (|2000| 
and Kotov, Churazov, & Gilfanov (2001) have suggested 
that steps in lag vs. frequency could be related to the tran- 
sition between two variability components, since the steps 
seem to occur in between the 'bumps' or broad Lorentzians 
in the PSD. If each component shows a different lag, one 
would expect a stepped lag vs. frequency dependence con- 
sisting of discrete lags with relatively weak lag vs. frequency 
dependence where a single component dominates the vari- 
ability, but with sharp steps occuring in the transition re- 
gion between components. This picture may fit with the 
SWIFT J1753. 5-0127 2009 data, but the signal-to-noise is 
such that it is difficult to find a precise match between the 
step in the lag and the transition between the low and high- 
frequency components in the PSD (which seems to occur 



5 CONCLUSIONS 

The soft X-ray coverage of the EPIC-pn camera onboard 
XMM- Newton has allowed us to carry out a comprehen- 
sive study of the spectral-timing properties of the BHXRB 
SWIFT J1753.5-0127 in the luminous 2009 hard state, and 
to perform a comparison with the fainter 2006 hard state. 

In our previous study of a sample of hard-state 
BHXRBs (Uttley et al. 2011), we showed that (i) soft (disc) 
X-ray variations lead hard (power-law) variations up to fre- 
quencies of about 1 Hz, and (ii) at the same low frequen- 
cies, the fractional variability in the disc is larger than 
the variability observed in the power-law component. Mass- 
accretion rate fiuctuations that propagate through the disc 
would produce variability on frequencies corresponding to 
the viscous time-scale where they are generated, and even- 
tually reach a centrally-concentrated source of power-law 
photons. This basic picture would both explain the hard-to- 
soft lags and give a reasonable explanation for the excess 
fractional variability seen in the disc as compared to the 
power-law. 

The 2009 observation shows a breadth of previously un- 
known observational features, which we summarise as fol- 
lows. 

Firstly, the variability amplitudes in the PSD are 
strongly dependent on energy, with the harder band varying 
more than the softer bands. The explanation for this effect is 
seen in the fractional covariance spectrum, where the bright 
disc black-body component that is seen in the spectrum now 
appears suppressed in fractional-variability terms with re- 
spect to the power-law, at all frequencies considered. This is 
opposite to what was found in less-luminous hard states in 
Uttley et al.| ( |2011p , where the disc was more variable than 
the power-law component. The higher luminosity (and thus, 
accretion rate) seen in this observation could imply that the 
disc emissivity in the observed bands is less centrally peaked 
than in less-luminous hard states. In this case, a fluctuation 
that is produced at an inner disc radius (e.g. within 1.5 nn) 
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will not modulate the emission from outer radii and this 
will effectively reduce the observed fractional variability of 
the disc as a whole. This is consistent with the shape of the 
energy-dependent lags. 

Secondly, wc have found that the low frequency (< 
0.5 Hz) disc variability below < 0.6 keV appears partly un- 
correlated with the rest of the emission produced by both 
the disc and power-law components. An explanation for this 
could exist if this incoherent accretion variability is pro- 
duced at larger radii than the coherent variability, affecting 
the softest emission that is more extended while being un- 
able to propagate to the innermost part of the disc and the 
power-law emitting region due to viscous damping. 

Finally, the last feature to highlight is an up-turn in 
the frequency-resolved time lags above ~ 0.5 Hz. This could 
be produced if, e.g. a relatively central secondary extended 
corona that sandwiches the disc - and thus reduces the ob- 
served hard-to-soft lags caused by disc propagation varia- 
tions - becomes swamped by a coherent source of variations 
associated with the high-frequency peak in the PSD (per- 
haps an even smaller unstable radius in the disc). 

By identifying a number of spectral-timing signatures 
with the accretion disc, this work suggests that a wealth 
of information regarding the properties of accretion flows 
can be discovered by extending spectral-timing measure- 
ments to softer X-rays and using techniques that allow bet- 
ter spectral-resolution. By combining these powerful tech- 
niques with more detailed modelling, a comprehensive un- 
derstanding of the variability may help to unravel many of 
the mysteries underlying accretion physics and associated 
phenomena such as the nature of state transitions and the 
production of the jets seen in the hard states of accreting 
black holes. 
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